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Abstract  of  Dissertation  Presented  to  the  Graduate  Council 
In  Partial  Fulfillment  of  the  Requirements  for  the 
Degree  of  Doctor  of  Philosophy 

MAGNETIC  FIELD  SOLUTIONS  FOR  ONE  CLASS 
OF  NONLINEAR  BOUNDARY  VALUE  PROBLEMS 

by 

Fred  Oliver  Simons,  Jr. 

June  22,  1965 

Chairman:  Dr.  W.  E.  Lear 

Major  Department:  Electrical  Engineering 

Analytical  d.c.  magnetic  field  solutions  to  nonlinear 
boundary  value  problems,  which  result  when  ferromagnetic 
materials  are  placed  In  magnetic  fields,  are  conspicuously 
missing  from  the  literature.  The  work  presented  in  this 
dissertation  presents  one  approach  which  yields  acceptable 
results  for  many  of  these  types  of  boundary  value  problems. 

Any  attempt  to  solve  the  general  anisotropic  non- 
linear boundary  value  problem  would  most  likely  be  futile. 
Therefore,  a system  of  analytical  techniques  was  applied  to 
a simplified  model,  which  is  a good  approximation  to  many 
different  types  of  the  original  more  difficult  problem.  The 
simplified  model  essentially  consisted  of  approximating  the 


lonhomogeneous  magnetic 


anisotropic  magnetic  medium  as  a n 
medium.  Although  the  nonhomogeneous  boundary  value  problem 
Is  still  nonlinear.  It  Is  shown  that  solutions  can  be 
obtained . 

The  B-H  loop  for  annealed  low  oarbon  steel  was  used  to 
obtain  the  nonhomogeneous  magnetic  characteristics.  The  only 
quantity  needed  from  the  B-H  loop  was  the  magnetic  permea- 
bility ( /£  ),  which  wsb  shown  to  be  representable  as  a func- 
tion of  magnetic  flux  density  (B). 

Maxwell's  equations  for  the  time -Invariant  case  were 
used  to  derive  the  two  equations  which  describe  the  magnetic 
fields  for  the  nonhomogeneous  problem.  These  equations  are 
put  into  a more  useful  form  by  defining  and  Introducing 
vector  magnetic  potential  Into  the  equations. 

In  order  to  Illustrate  that  the  equations  are  valid,  a 
particular  two-dimensional  boundary  value  problem  Is  posed 
and  solved.  Three  major  steps  were  required  to  obtain  the 
solutions.  First,  the  equations  for  vector  magnetic  potential 
had  to  be  put  Into  a form  which  applied  to  the  given  two- 
dimensional  problem.  The  magnetic  characteristics  had  to 
be  put  into  an  analytical  form  at  this  point  also.  Second, 


boundary  oonditiona  which  applied  to  the  problem  being 
treated  were  formulated  and  presented  in  a composite  form. 
Third,  finite  difference  schemes  were  devised  for  the 
equations  which  describe  the  problem.  These  schemes  were 
used  to  convert  the  problem  into  a computer  program  which 
was  solved  utilizing  the  IBM  709  computer. 

The  magnetic  field  computer  solutions  for  the  particular 
boundary  value  problem  were  interpreted  and  shown  to  yield 
the  expected  results.  It  was  easy  to  see  that  as  the  sample 
material  was  exposed  to  continuously  stronger  fields,  it 
became  saturated. 

Thus,  it  can  be  concluded  that  the  major  contribution 
of  this  dissertation  is  that  new  techniques  have  been 
generated  which  can  be  applied  to  a whole  new  class  of  non- 
linear boundary  value  problems.  This  class  can  be  described 
as  those  d.c.  magnetic  nonlinear  boundary  value  problems 
which  result  when  the  magnetic  medium  is  either  nonhomogeneous 


approximately  nonhomoger 


CHAPTER 


INTRODUCTION 

In  recent  years  microwave  technology  has  been  Involved 
In  a continual  state  of  growth  and  advancement  due  to  the 
demands  created  by  many  phaseB  of  communication.  In  the 
process  of  meeting  these  new  demands,  many  related  problems 
are  continually  confronted. 

One  such  problem,  which  was  brought  to  the  attention 
of  the  author  by  Dr.  A.  D.  Sutherland  (of  Sperry  Microwave 
Tube  Division  of  the  Sperry  Rand  Corporation  of  Oalnesvllle, 
Florida),  Is  that  of  determining  the  magnetic  field  dis- 
tribution In  the  electron  gun  region  of  traveling  wave  tubes. 
High  frequency  and  high  power  specifications  require  highly 
convergent  electron  beams . This  condition  results  since  high 
frequency  dictates  smaller  physical  dimensions  of  the  beam 
while  high  power  would  require  higher  beam  energy  densities. 
If  the  beam  trajectories  coincide  with  the  magnetic  flux 
lines,  the  electron  gun  problem  is  reduced  to  an  electro- 
static electron  gun  design.  Since  the  d.c.  magne 


stic  focusing 


circuit  is  terminated  in  the  electron  gun  neighborhood,  the 
resulting  magnetic  field  has  the  required  high  degree  of 
convergence.  Thus,  if  the  magnetic  field  distribution  is 
known  and  if  the  beam  trajectories  are  made  to  coincide 
with  the  magnetic  field  lines,  the  electron  gun  design 
is  greatly  simplified.  Heretofore,  the  magnetic  field  in 
the  gun  region  has  been  found  from  iron  filing  patterns. 
However,  this  method  not  only  requires  a previously  con- 
structed focusing  system,  but  it  also  becomes  increasingly 
unreliable  for  the  smaller  dimensions  that  are  being  en- 
countered. In  view  of  these  points,  it  would  be  desirable 
if  an  analytical  method  of  analysis  could  be  derived. 

Thus  the  problem  posed  is : 

Given  the  geometry  of  a magnetic  circuit 
with  a magnetic  focusing  material,  which  13 
not  only  nonhomogeneous  but  also  anisotropic, 
to  find  the  magnetic  field  configuration  to 
a sufficient  degree  of  accuracy  by  some 

It  should  be  noted  that,  although  the  problem  as 
posed  above  is  an  outgrowth  of  microwave  advancements,  any 
solution  techniques  generated  will  apply  to  a much  more 
general  class  of  nonlinear  boundary  value  problems.  In 


general,  the  d.e.  magnetic  field  problem  is  a nonlinear  and 
nonhomogeneous  boundary  value  problem  with  an  anisotropic 
magnetic  medium  present. 

Before  a solution  is  attempted,  the  anisotropic 
magnetic  properties  will  be  approximated  as  nonhomogeneous 
properties.  Arguments  are  presented  which  validate  the 
approximation.  The  partial  differential  equations  which 
apply  to  each  of  the  media  involved  (free  space,  a current 
region,  and  a nonhomogeneous  magnetic  medium)  are  then 
derived  from  Maxwell's  equations.  Next,  a two-dimensional 
mesh  is  presented  which  applies  to  a particular  boundary 
value  problem.  Finite  difference  schemes  are  devised  for 
the  partial  differential  equations  and  applied  to  the  chosen 
boundary  value  problem.  Finally,  solutions  are  obtained  by 
programing  the  problem  on  the  IBM  709  high-speed  digital 
computer  and  the  results  are  interpreted  to  demonstrate  that 
the  new  solution  techniques  are  valid. 


MAGNETIC  CHARACTERISTICS 


In  the  most  general  sense.  It  would  be  virtually 
Impossible  to  depict  the  complete  anisotropic  nature  of  a 
ferromagnetic  material.  The  behavior  of  a point  in  the 
material  is  not  only  dependent  on  the  value  of  the  applied 
field,  but  it  also  depends  on  the  previous  and  the  present 
states  of  all  the  other  points  in  the  material.  The  most 
fundamental  magnetic  characteristic  is  a simple  hysteresis 
loop,  which  is  more  commonly  referred  to  as  a B-H  loop. 

The  B-H  loop  for  a particular  sample  of  annealed  low  carbon 
steel  is  Shown  in  Figure  1 (l) .*  The  characteristics  of 
this  sample  will  be  used  for  the  particular  boundary  value 
problem  that  will  be  treated. 

Notice  that  the  state  of  any  differential  volume 
within  the  sample  could  be  represented  by  any  point  on  the 
boundary  or  inside  the  B-H  loop.  A very  simple  approximation 


•Numbers  in  parentheses  denote  entries  in  the  List  of 
References . 


H (oersteds) 

Figure  1.  Magnetic  Characteristics  for  Annealed 
Low  Carbon  Steel 


to  the  magnetic  characteristics  would  be  a single  valued 


function  of  magnetic  flux  density  vs.  magnetic  field 
intensity.  Such  a characteristic  is  shown  as  a dotted  line 
in  Figure  1 . 

There  are  two  reasons  why  this  type  of  characteristic 
would  be  a good  approximation.  First,  the  B-H  loop  is 
narrow  enough  so  that  large  errors  do  not  result,  and  these 
errors  are  seen  to  be  much  smaller  when  the  sample  is  driven 
near  saturation.  For  example,  errors  of  the  approximation 
are  less  than  ten  percent  for  any  point  in  the  sample 
where  B is  greater  than  8 kllogauss . In  magnetic  focusing 
structures,  materials  are  usually  operated  in  or  near 
saturation.  Second,  if  suoh  a simplified  characteristic 
were  not  valid  the  magnetic  field  would  not  repeat  when 
the  magnetic  circuit  was  energized  at  different  times  and. 
to  varying  degrees.  This  behavior  has  not  been  observed 
in  practice,  or  if  it  has,  it  has  been  negligible. 

In  summary,  the  anlstroplc  ferromagnetic  material 
characteristics  will  be  approximated  as  nonhomogeneous 
magnetic  material  characteristics.  This  approach  is  not 
only  useful  for  the  above  reasons,  but  it  will  be  shown 


smatically  expedient. 


the  magnetic  medium  is  usually  energized  and 
deenergized  in  one  direction  only,  the  single  valued 
characteristic  that  is  used  is  assumed  to  lie  near  the  top 
of  the  B-H  loop.  The  curve  was  also  chosen  from  the  stand- 
point of  getting  a good  fit  while  at  the  same  time  having  a 
form  easy  to  treat  mathematically.  The  curve  chosen  has  the 
form  given  by 

B = 4H  0 = H 5 1 (2-1) 

B = 10  - 6e"  ^3^-  1 5 H ' 13.45  (2-2) 

B = 9.998  + 0.001  H H 2 13.45  (2-3) 

where  B is  kilogauss  and  H is  in  oersteds.* 

The  magnetic  permeability,  which  is  defined  by 

B/H  (2-4) 

and  which  is  derived  from  equations  (2-1)  through  (2-3),  is 
given  by 

H - 16  7T  (10‘4)  0 2 B — 0.4  (2-5) 

H „ 40  7TB  (IQ'4) 

1 + i.5  in  (2c|)  0.4  2 B — 0.999&  (2-6) 

LL  = 4 7T  B (10~7) 

(B  - 0.9998)  + 13.45( .0001) 

B - 0.9998  ( 2-7) 

•Unless  otherwise  specified,  all  symbols,  which  are 
defined  in  the  List  of  Symbols,  are  in  the  MKS  syste 
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The  second  breakpoint  for  B was  chosen  such  that  the 
slope  of  the  B-H  characteristic  was  Mo.  This  was  done 
since  just  above  saturation,  the  B-H  characteristic  becomes 
a straight  line  with  slope  Mo- 
lt Is  shown  later  that  the  derivative  of  V ( l>~  1/fA) 
with  respect  to  B (£/')  Is  needed  In  the  problem  formulation. 
Therefore,  It  is  presented  at  this  time. 

0 0 5 B 2 0.4  (2-8) 


V'- 


1.5  - 1_+  1.5  In  ut§ 

bTi-bT  b2  52 


V'J-# 5 do5) 


0.4  2 B 2 0.9998  (2-9) 


B 2 O.9998  (2-10) 


In  summary  the  necessary  magnetic  characteristics 
of  a particular  ferromagnetic  material  have  been  formulated 
and  presented  in  equations  (2-6)  through  (2-10).  The 
functions  described  by  these  equations  are  plotted  in 
Figure  2.  The  digital  computer  program  that  was  used  to 
obtain  the  second  breakpoint  of  B and  to  solve  for  /J.  and  V 
as  functions  of  B is  given  in  Appendix  A. 


B (weber) 

Figure  2.  Plots  of  fi  and  U'  Normalised  with 
Respect  to  their  Maximum  Values 


CHAPTER  3 


MAONETIC  FIELD  EQUATIONS 

There  are  certain  restrictions  Tor  the  particular 
boundary  value  problem  that  will  be  treated.  The  first 
of  these  Is  that  the  problem  Is  time  Invariant;  therefore. 
Maxwell's  equations  are  decoupled  and  it  is  only  necessary 
to  deal  with  the  two  magnetic  field  relations.  Further, 
the  problem  will  be  for  an  Isotropic  nonhomogeneous  magnetic 
medium  In  the  presence  of  free  space  and  current  carrying 
conductors.  Of  course  In  free  space  and  In  the  more  common 
nonferromagnetic  conductors, 

u-  Ho  ■ (3-D 

In  the  nonhomogeneous  magnetic  medium,  p.  and  \j‘ 
will  be  treated  as  Implicit  functions  of  space  coordinates. 
With  these  restrictions,  the  two  magnetic  field  relations  of 
Maxwell's  equations  become 


V-  5 - o 

yx  H = T 


10 


(3-2) 

(3-3) 


auxiliary  relatlc 


B - M B (3-4) 

is  understood  to  be  defined  for  M a function  of  space 
coordinates.  Equation  (3-4)  can  be  used  to  change  the 
form  of  (3-3)  to  get  the  set  of  equations  (3-5a)  and 
(3-5b). 

V-  E = 0 (3-58) 

VX  ( V 5)  - T (3-5b) 


V 


(3-6) 


Equations  (3-5a)  and  (3-5b)  agree  with  the  magnetic  field 
equations  derived  by  Smythe  (2)  for  the  same  conditions. 

The  equations  (3-5a)  and  (3-5b)  could  be  applied  in 
their  present  form;  however,  it  can  be  shown  that  a new 
vector  quantity,  the  vector  magnetic  potential  ff,  is  easle 
to  utilize.  The  defining  equations  for  vector  magnetic 
ntial  are  given  by  (3,4); 


12 


V *-°  (3-7a) 

Vxff-  B . (3.7b) 

When  equation  (3-7b)  la  aubBtltuted  Into  equations  (3-5a) 
and  (3-5b),  the  relations  obtained  are 

V-  V X - 0 (3-8a) 

Vx(yyxX)-T.  (3-8b) 

If  solutions  for  A can  be  obtained  from  equations  (3-8a)  and 
(3-8b),  then  equation  (3-7b)  can  be  used  to  find  the  magnetic 
field.  Note  that  equation  (3-8a)  is  automotlcally  satisfied 
by  a vector  identity  (5).  Thus,  It  is  only  necessary  to 
solve  the  equation  (3-8b)  for  A.  This  simplification  In 
the  equations  was  obtained  by  using  the  defining  relations 
(3-7a)  and  (3-7b). 

The  solution  for  the  magnetic  field  boundary  value 
problem  that  has  been  outlined  has  been  reduced  to  that  of 
finding  solutions  to  equation  (3-8b),  which  can  be  expanded 
In  the  form  of  equation  (3-9)  or  (3-10)  by  using  another 


Identity  (5). 


2/v*  v**+  vyx(VxA)  = T 

Vx  yx  a = /ir  - /iVy  X ( Vx  a)  . 


(3-9) 


(3-10) 


If  one  other  vector  Identity  is  applied  to  equation  (3-10) 
(5),  then  the  final  relation 


is  obtained.  Equation  (3-11)  is  a general  relation  and  as 
such  is  applied  to  free  space,  magnetic  medium,  and  current 
regions . 

The  equation  for  vector  magnetic  potential  in  free  space 
is  given  by 


and  it  is  obtained  from  equation  (3-11)  by  recognizing  that 


V2  A . /iVyx(VXA)  - /J.T 


(3-11) 


V*  A - 0 


(3-12) 


1=0 


(3-13) 


n-  /J-o  or  Vy  - o 


(3-14) 


The  equation  for  the  vector  magnetic  potential  in  a 


mferromagne 


conductor 


I is  given  by 
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V2  A - - Mo  T , (3-15) 

Vy  - 0 . (3-16) 

The  equation  for  the  vector  magnetic  potential  In  the 
magnetic  focusing  material  Is 

V2  *-  MVyX(VXA)  (3-17) 

and  is  obtained  from  (3-11)  slnoe  I - 0. 

In  summary,  three  equations 

V2  A = 0 • (3-18) 

V2A  --MT  (3-19) 

V2A-  ^VyX(VXA)  (3-20) 


are  derived  which  describe  the  vector  magnetic  potential 
In  free  space.  In  a current  carrying  non-magnetic  conductor, 
and  In  a nonhomogeneous  magnetic  medium  respectively.  These 
relations  will  be  used,  with  formulated  finite  difference 
schemes.  In  a computer  program  to  solve  a particular  non- 
linear bo 


nundary 


CHAPTER 


THE  SOLUTIONS  TO  A PARTICULAR 
BOUNDARY  VALUE  PROBLEM 

The  ba3lc  relations  which  are  needed  to  solve  one 
class  of  nonlinear  magnetic  field  boundary  value  problems 
have  been  derived  and  presented  in  Chapters  2 and  3.  In 
the  present  chapter  these  techniques  will  be  applied  to  a 
particular  problem  which  will  illustrate  that  the  derived 
relations  are  correct.  It  should  be  noted  that  the  appli- 
cation of  these  relations  constitute  a major  portion  of 
the  difficulty  in  solving  problems  of  the  nature  outlined. 

A problem,  which  is  simple  enough  to  solve  without 
undue  expense  and  time  and  at  the  same  time  complex  enough 
to  illustrate  nonlinear  solutions,  is  presented  in  Figure  3. 
It  should  be  pointed  out  that  this  problem  was  the  last  of 
three  that  were  seriously  considered.  As  simple  as  it  is, 
the  problem  required  one  and  one-half  hours  computer  time  on 
the  IBM  709  computer.  A simple  description  of  the  problem 
will  serve  to  clarify  Figure  3- 

15 
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Figure  3. 


The  Cross  Section  Geometry  of  the  Non- 
linear Boundary  Value  Problem  with  a 
Uniform  Gap  Field  E0  Applied  Between 
Infinite  jX  Pole  Faces 


s\\\\\\\\\\>X\\\\\\\^\\\\ 
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What  Is  the  two  dimensional  magnetic  field 
configuration  produced  when  a bar  of  square 
cross  section  with  nonhomogeneous  magnetic 
characteristics  Is  Inserted  in  a uniform 
gap  field? 

The  boundary  conditions  for  the  problem  depicted  In  Figure  3 
will  be  given  In  Section  4-1.  The  equations  that  apply  to 
this  problem  will  be  presented  In  Section  4-2.  Finally,  the 
results  obtained  will  be  presented  In  Section  4-3. 

Section  4-1  — Boundary  Conditions 

The  boundary  conditions  could  be  formulated  for  the 
problem  as  It  Is  presented  In  Figure  3.  However,  due  to 
symmetry.  It  will  only  be  necessary  to  treat  one  quadrant  of 
the  geometry.  Thl3  quadrant  is  presented  In  Figure  4 along 
with  one  typical  mesh  point. 

In  order  to  formulate  the  boundary  conditions,  the 
integral  relation  (6) 


should  be  considered.  Note  that 


1 - V*  H 


(4-2) 
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Boundary  Value  Problem  with  a Typical  Mesh 
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from  one  of  Maxwell's  equations,  equation  (3-3).  Since  H 
has  only  x and  y components  and  13  only  a function  of 
x and  y,  A has  only  a z component  (A  = Aaz)  and  Is  a function 
of  x and  y.  Also,  since  the  x-z  plane  Is  a constant  flux 
plane,  then  it  can  be  chosen  as  a zero  potential  reference 
plane.  Further,  from  the  symmetry  chosen,  the  boundary 
condition  for  A at  the  y-z  plane  becomes 


Finally,  the  boundary  condition  for  the  outer  extremity  (X) 
of  the  problem  can  be  obtained  from  the  function 


since  A Is  only  a function  of  y at  x = X.  Note  that  the 
field  is  uniform  at  X;  then 


which  is  the  final  boundary  condition  needed.  This  follows 


$=A 


(H-5) 


B0Y 


for  y = Y 


(4-6) 


which  may  be  obtained  simply  from  equation  (4-5). 

For  a composite  reference,  the  previously  derived 
boundary  conditions  are  complied  and  presented  here. 


A = 0 y = 0 

HL  0 v n 

A = B0y  X = X 

A = B0Y  y = Y 


0 2 x - X (4-7a) 

o — y — Y (4— 7l>) 
o iy  2 Y (4-7c) 
0 ' X ; X (4-7d) 


Section  4-2  --  Equation  Formulation 

The  equations  that  are  used  to  solve  the  boundary  value 
problem  will  first  be  presented  In  the  partial  differential 
equation  form  and  then  In  the  finite  difference  form. 

Partial  differential  equations  can  be  obtained  by  a 
direct  application  of  the  problem  restrictions  to  equations 
(3-18),  (3-19) ( and  (3-20).  Equation  (3-20)  will  be  treated 
first  and  then  the  applicable  equations  may  be  compiled  by 
Inspection.  Since  A has  only  a z component  and  Is  only 


function  of  x and  y,  then 


(4-8) 
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Vs  A. 


[a2 a 
[o* 


a2fl 


'K 


Due  to  the  fact  that  /I  and  V for  the  time  being  are  con- 
sidered to  be  Implicit  functions  of  space  coordinates, 

(x  and  y)  then 


fe-Sr-oF-Sth-'*-9’ 

At  this  point  the  magnetic  medium  characteristics  can 
be  taken  Into  aooount.  The  magnetic  flux  density 

B - V X A (4-10) 

has  a magnitude  given  by 

b ■ 

The  quantities  /i  and  V were  shown  to  be  representable 
as  functions  of  B alone  by  equations  (2-5)  through  (2-10). 


Thus  these  relations  fc 
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Qv 

(4-12a) 

Ov 

■ 

(4-12b ) 

The  partial 

derivatives  of  B with  respect  to 

x and  y are 

obtained  frc 

>m  equation  (4-11) 

<3b 

,i  [a*  0!a  + ^a  dh_] 
B [37  Qy  axayJ 

(4-13a) 

and 

-Qs- 

ay 

I [d-A  Ogft  X (9a  a2A 

B lax  axay  ay  "a^j 

(4-13b) 

These  latter  relations  Illustrate  that  U and  V satisfy 
the  condition  specified  earlier;  l.e.,  they  are  representable 
as  Implicit  functions  of  space  coordinates. 

Equations  (4-8)  through  (4-13)  are  the  relations 
ed  to  simplify  equation  (3-80)  to 
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dh  + 52a  _ . uv'  [ _^a  I agA 

"sf  » [ a- 1 a- 

. a>  aa»  ' ..  a>  I a»  agA 
a'S^a7)  a’la-d'a' 

* a.  a »«\1 
3»c5>sJJ  - 


Another  form  of  (4-14)  Is 

S4*^'(SC] 

(Si)1 

+ zA/a*  ^a  agA  = o . 

b ay  axay  (*-»> 


This  relation  serves  better  to  illustrate  the  nonlinear 
nature  of  the  problem  since  products  and  squares  are 
displayed  together. 

With  equations  (4-8)  and  (4-14),  the  equations  which 
are  necessary  to  describe  the  particular  boundary  value 


given  by 


(1.-168) 
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„ < 3 a a^A  1 

a»  a*  ey  / 


where  equation  (4-l6a)  refers  to  the  free  space  region 
while  (4-l6b)  refers  to  the  magnetic  region.  If  there 
were  currents  present  in  the  problem  to  be  treated,  ther 
a third  relation  obtained  from  equation  (3-19)  would  be 


(.-IT) 


where  the  z component  of  B is  still  zero.  Although  this 
relation  is  not  needed,  it  is  worth  noting  that  equations 
(4-l6b)  and  (4-17)  can  be  used  to  identify  an  equivalent  z 


jomponent  of 
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(4-18) 


This  iterates  the  fact  that  magnetic  polarization  effects 
can  be  replaced  by  equivalent  currents  (6). 

The  boundaries  between  the  magnetic  medium  and  free 
space  constitute  an  abrupt  change  in  magnetic  permeability. 
The  normal  components  of  magnetic  flux  density  and  the 
tangental  components  of  magnetic  field  intensity  must  be 
matched  at  these  boundaries.  The  pertinent  equations  are 


Bni  - Bn2 
Htl  - Hta  . 

From  the  definition  of  5,  that  is 
VXA-B  , 

equations  (4-19)  and  (4-20)  can  be  transformed  to 


[0A  1 

[aA  I 

Lo^J, 

(4-19) 

(4-20) 

(4-21) 


(4-22) 
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The  first  of  these  two  relations  is  automatically  satis- 
fied since  A is  a continuous  function  throughout  the  mesh. 
The  second  equation  is  satisfied  when  the  first  order 
derivative  approximations  are  matched  on  each  side  of  the 
boundary.  The  value  of  A on  the  boundary  can  be  shown  to 
be  a function  of  adjacent  values  of  A that  are  not  on  the 
boundary  Itself.  The  finite  difference  scheme  applied  to 
equation  (4-23)  will  be  presented  along  with  other  finite 
difference  forms. 

All  the  equations  which  are  necessary  to  describe 
the  particular  boundary  value  problem  have  been  formulated . 
These  equations  (4-l6a),  (4-l6b)  and  (4-23)  can  be  put  into 
similar  finite  difference  forms  by  using  the  first  order 
approximations  (7).  These  are 


Wl,n  ~ Am-1 


(4-24a) 
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Am,n+1  - Am,n-1 

2 3 


(4-24b) 


[wvu 


oa».  > . 

~5*  * Vi 

dgA  - i [ . 

'9*3*  *"5*  'Wl-n+1  + 


(4-24o) 


(4-24d) 


(4-24e) 


When  equations  (4-24)  are  applied  to  equations  (4-l6a)  and 
(4-l6b),  the  results  are 


Vn 


, Arotl,n  * *»-l,n 


Am,n+1  *a,n-l 


(4 -25a) 


v»-  [»(i*-4£aj]'1  [ w 
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Am-l,n)  (*»fl,n  " Am.ljn)2  + 2 


i,n-l)  (Am,iHl  " Am,n-1^ 


+ (Ann-l,n  ■ Am-X,n^  (Am,r»l  ' Am,n-l)  (Amtl,nfl 


+ Ni-l.n-l  " Am-l,n+l  * Ann-l,n-l) 


)] 


(|t-25b) 


The  application  of  equations  (4-24)  to  equation  (4-23)  for 
a boundary  which  is  parallel  to  the  x-axls  would  yield  the 
relation 


A „ //m,n+l  Am,n-1  + Mm.n-1  ^.w-l 

m,n  A4m,n+1  + H-m.n-1  (4-25c) 

where  other  boundary  orientations  would  yield  similar 
relations . 

The  boundary  conditions  In  equation  set  (4-7)  and 
finite  difference  relations  in  (4-25)  constitute  the  necessary 
relations  that  are  needed  to  solve  the  problem  of  Figure  4. 
These  relations  were  converted  to  a digital  computer  program 
which  Is  presented  in  Appendix  B.  In  essence,  the  program 
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performs  the  function  of  applying  the  correct  finite 
difference  equation  to  each  mesh  point  In  an  Iterative 
fashion.  The  Initial  values  of  A,  before  relaxation, 
were  chosen  to  correspond  to  the  uniform  magnetic  field 


Section  4-3  — Results 

The  boundary  value  problem  Illustrated  In  Figure  3 
was  solved  four  different  times  with  the  computer  program 
in  Appendix  B.  The  four  oase3  are  enumerated  below. 

1.  The  boundary  value  problem  was  first  solved 

by  treating  the  magnetic  medium  as  homogeneous; 
i.e.,  IX  = 4000  and  V"  = 0 for  all  values  of 
B.  This  type  problem  Is  solved  In  the  liter- 
ature (7).  However,  this  solution  not  only 
serves  as  a check  on  the  equations,  but  It 
should  be  compared  to  the  computer  solutions 
for  other  cases.  It  should  be  noted  that  the 
computer  program  In  Appendix  B for  this  case 
had  to  be  modified  to  obtain  this  particular 
solution. 

2.  With  the  computer  program  of  Appendix  B,  a 
solution  was  obtained  with  a gap  field  of 

2 kllogauss  applied.  The  result  was  observed 
to  give  a maximum  field  of  3.94  kllogauss. 

Thus  the  nonlinear  effects  should  not  be  In 
the  solution  since  V'  - 0 for  B less  than 
4 kllogauss.  This  would  indicate  that  the 
results  for  cases  1 and  2 should  be  the 
same.  They  were  identical!  Why  then,  was 
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the  second  case  treated?  Because  it  illus- 
trates that  the  nonlinear  program  would 
revert  to  the  linear  solutions  when  the  iron 
is  not  driven  into  the  nonlinear  region. 

3.  The  third  case  was  for  a gap  field  of  8 
kilogauss.  This  corresponds  to  a high 
degree  of  nonlinearity. 

4.  The  fourth  and  final  case  was  solved  with 
a gap  field  of  10  kilogauss,  which  is 
approximately  the  saturation  flux  density 
for  the  ferromagnetic  sample  treated . 

In  order  to  illustrate  these  cases,  the  results  obtained 
for  cases  1 and  2 are  presented  in  Figure  5 while  the  results 
for  cases  3 and  4 are  presented  in  Figures  6 and  ^ respec- 
tively. Each  of  these  magnetic  field  plots  was  obtained 
when  constant  contours  of  the  z component  of  vector  magnetic 
potential  were  sketched  directly  on  the  two  dimensional 
arrays  of  numbers  provided  by  the  computer  output.  Then 
these  plots  were  transferred  to  Figures  5,  6,  and  7. 

There  are  primarily  two  observations  that  serve  to 
validate  the  techniques  for  solving  the  particular  nonlinear 
problem  that  was  treated.  First,  as  the  magnetic  field 
increased,  the  magnetic  flux  lines  began  to  distort  within 
the  magnetic  material.  This  indicated  that  M was  becoming 
nonhomogeneous  as  the  magnetic  material  began  to  saturate. 
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Figure  5. 


Magnetic  Flux  Line  Solutions  for 
the  Nonlinear  Boundary  Value  Problem 
with  a Gap  Field  of  2 Kllogause 


32 


Figure  6. 


Magnetic  Flux  Line  Solutions  for 
the  Nonlinear  Boundary  Value 
Problem  with  a Gap  Field  of 
8 Kllogauss 
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Figure  7. 


Magnetic  Flux  Line  Solutions  for 
the  Nonlinear  Boundary  Value 
Problem  with  a Gap  Field  of  10 
Kllogau33 


especially 


The  second  point,  which  adds 


credenoe  to  the  solutions  obtained,  is  the  ratio  of  the 
field  Just  inside  the  magnetic  material  to  that  Just  outside 
the  material . This  ratio  was  obtained  at  the  plane  of 
symmetry  where  S is  parallel  to  the  boundary.  The  resulting 
ratio  parameters  have  been  computed  and  enumerated  below 
for  reference. 

1.  The  ratio  for  cases  1 and  2 is  given 

by  R - 10  . (4-26a) 

2.  The  ratio  for  ease  3 is  given  by 

R - 3.75  (4-26b) 

where  the  gap  field  was  8 kilogauss. 

3.  The  ratio  for  case  4 is  given  by 

R = 3.00  (4-26c) 


where  the  gap  field  was  10  kilogauss. 

These  results  are  exactly  as  expected;  i.e.,  a3  the  gap 
field  is  Increased,  the  magnetic  material  shows  definite 
trends  toward  saturation. 


CHAPTER  5 


CONCLUSIONS 

Before  any  conclusions  are  presented,  it  would  be 
appropriate  to  summarize  the  work  that  has  been  presented. 

The  original  problem  posed  wa3,  "can  solutions  be  obtained 
for  the  general  anisotropic  nonlinear  boundary  value 
problem?"  The  form  of  the  magnetic  characteristics 
allowed  the  anisotropic  magnetic  medium  to  be  approximated 
as  a nonhomogeneous  medium;  i.e.,  fJ.  was  representable 
as  a function  of  space  coordinates.  After  the  vector 
magnetic  potential  was  Introduced,  it  was  applied  to 
Maxwell's  d.c.  magnetic  equations  to  obtain  a relation 
which  was  applicable  to  a current  carrying  nonhomogeneous 
magnetic  medium.  This  general  relation  and  the  magnetic 
medium  characteristics  were  used  to  obtain  the  equations 
that  would  apply  to  free  space,  current  regions,  and  non- 
homogeneous magnetic  regions.  Then  boundary  conditions 
were  specified  for  a particular  problem,  and  this  problem 
was  programed  and  solved  on  the  IBM  709  computer.  The  results 
obtained,  which  were  positive,  were  presented  In  Section  4-3. 
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The  first  and  most  Important  oonoluslon  to  be  drawn 
Is  that  the  new  techniques  applied  to  this  class  of  non- 
linear boundary  value  problems  are  valid.  This  conclusion 
wa3  established  by  first  treating  a linear  problem  which  was 
simple  enough  to  give  known  solutions,  and  then  determining 
the  effects  of  driving  the  magnetic  material  Into  the 
nonlinear  region.  Hie  saturation  effects  observed  were 
exactly  those  expected. 

The  second  conclusion  concerns  the  complexity  of  the 
problem.  Although  the  problem  was  treated  by  approximating 
the  anisotropic  medium  as  a nonhomogeneous  magnetic  medium, 
this  does  represent  a step  in  the  direction  of  treating  a 
whole  new  class  of  nonlinear  boundary  value  problems.  In 
past  work  only  those  problems  Involving  homogeneous  magnetic 
media  have  been  solved.  Of  course,  the  usefulness  of  the 
techniques  will  be  determined  by  how  well  the  nonhomogeneous 
model  approximates  the  anisotropic  problem. 

There  are  certain  reasons  why  such  a simple  problem 
was  treated.  First,  a problem  was  chosen  for  the  sole 
purpose  of  proving  or  disproving  the  techniques  applied  to 
this  class  of  nonlinear  boundary  value  problems.  Naturally, 
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it  was  desirable  to  use  a problem  which  was  easy  to  interpret. 
Second,  as  simple  as  this  problem  was  in  form,  the  applica- 
tions of  the  techniques  were  quite  involved.  Approximately 
two  hours  of  computer  time  were  needed  to  complete  the 
program  and  solutions. 

The  techniques  that  have  been  presented  are  probably 
applicable  to  the  originally  posed  question,  which  was  to 
find  the  fringing  magnetic  field  in  the  gun  region  of 
traveling  wave  tubes.  This  conclusion  is  supported  by  the 
fact  that  the  focusing  material  of  a traveling  wave  tube 
has  a narrow  B-H  loop  as  in  Figure  1 . Also  the  field  plots 
needed  are  not  in  the  immediate  neighborhood  of  the.  iron 
which  is  where  the  errors  would  be  largest. 

A final  conclusion  that  can  be  drawn  is  that  the 
techniques  that  are  presented  here  warrant  a detailed 
theoretical  and  experimental  study.  Two  reasons  can  be 
presented  which  verify  this  conclusion.  First,  with  all 
the  work  and  attention  that  has  been  focused  on  the  homo- 
geneous problem,  it  is  certainly  Justifiable  to  treat  in 
detail  the  more  descriptive  nonhomogeneous  problem.  Second, 
an  extensive  study  is  needed  to  extend  the  techniques  to 
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axially  symmetric  problems,  three-dimensional  problems,  and 
other  more  general  cases. 

As  a closing  remark,  It  can  be  concluded  that  the 
major  contribution  of  this  dissertation  Is  that  new  techniques 
have  been  generated  which  can  be  applied  to  a whole  new  class 
of  nonlinear  boundary  value  problems.  This  class  can  be 
described  as  those  d.c.  magnetic  nonlinear  boundary  value 

nonhomogeneous  or  approximately  nonhomogeneous . 
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APPENDIX  A 


The  computer  program  which  Is  presented  here  was  used 
to  solve  for  M and  V as  a function  of  B.  The  defined 
symbols  resemble  the  quantities  they  represent  closely 
enough  to  allow  Interpretation  of  the  program  without  a 
detailed  cross-reference. 
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APPENDIX  B 


The  computer  solution,  which  was  used  to  solve  the 
particular  boundry  value  problem  in  Chapter  it,  is  presented 
at  this  time  for  reference.  A complete  explanation  of  the 
manner  in  which  the  problem  was  programed  will  not  be 
presented  for  the  salce  of  brevity,  but  those  Interested  in 
the  program  in  detail  would  not  have  any  difficulty  in 
interpreting  the  program  in  its  present  form. 
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